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1. Introduction 

The response of a Brownian particle to an external force is one of the paradigms of 
statistical mechanics. Brownian motion is usually described as a Wiener process and 
a full mathematical analysis is possible based on the solution of the corresponding 
Fokker-Planck equation. However, a genuine Gaussian white noise does not exist 
in real systems, and the importance of other types of noises has long been clear. 
Such noises have been studied in great detail in zero-dimensional systems, and their 



specific properties are known to have a profound influence on the behavior of these 
systems The effect of the color of the noise on noise- induced transitions 

continues to be documented 0E1E]) and has been found to be quite dramatic since 
it can alter the type of transition and lead to new re-entrance phenomena. Other 
noise-induced effects have also been found to be sensitive to the correlation time 
of the noise, e.g., stochastic resonance synchronization [7], and transport in 
Brownian ratchets 

The two most commonly discussed examples of colored noise are the Ornstein- 
Uhlenbeck process and the dichotomous Markov process [H]- The latter has the 
great advantage that analytic results can often be derived. Nevertheless an essential 
technical difficulty restricts almost all the results obtained in the past to dynamics 
with no fixed points, or exclusively with stable fixed points, see for example [2]. 
We have recently made progress toward overcoming the technical difficulty that 
appears in the presence of unstable fixed points by identifying the source of spurious 
divergences that arise in the usual analytic approaches to the problem 9 , lOj, and 
are now in the position to consider these cases as well. 

The purpose of this paper is to apply the general analytic procedure we de- 
veloped in Ref. jlOj to two different stochastic phenomena driven by dichotomous 
Markov noise, namely Stokes' drift |TJ and hypersensitive transport [T2"lll3j . Sec- 
tion [21 deals with Stokes' drift as a rocking ratchet problem and Sec. |3| with hyper- 
sensitive transport. A short summary is presented in Sec. Some of the lengthy 
calculations have been collected in the Appendices. All the analytical results are 
supported by numerical simulations of the corresponding dynamics of an ensemble 
of 20000 particles, with random initial positions, sampled at 100 different times in 
order to build up the histogram of the stationary probability density and, from it, 
to compute the value of the asymptotic mean velocity. 

2. Stochastic Stokes' Drift as A Rocking Ratchet 

A longitudinal wave traveling through a viscous fluid imparts a net drift to the 
suspended particles, an effect known as Stokes' drift. The deterministic effect (that 
does not account for the stochastic fluctuations or perturbations in the system) 
has been extensively studied in various practical contexts ranging from the motion 
of tracers in meteorology and oceanography .14 to doping impurities in crystal 
growth |15|. The deterministic drift has a simple intuitive explanation, namely that 
the suspended particle spends a longer time in the regions of the wave-train where 
the force due to the wave acts in the direction of the propagation of the wave than 
the time it spends in those regions where the force acts in the opposite direction. 
Therefore, the particle is driven on average in the direction of wave propagation. As 
a simple clarifying example , consider the dynamics of an overdamped particle 
forced by a traveling square wave with velocity v and wavelength L, 

x = f{x-vt), (1) 

where x(t) is the position of the overdamped particle at time t, and / is the periodic 
forcing due to the longitudinal wave traveling at the speed v and of wavelength L. 
Suppose that the particle is entrained with a force f = bv when in a crest part of 
the wave, and / = — bv when in a trough part. The time spent in a crest part, 



L/[2(l — b)v], is larger than in the though part, L/[2(l + b)v), resulting in a net 
deterministic drift velocity vq = b 2 v of the particle. 

Recent studies ^]E3E] show the importance of stochastic effects on Stokes' 
drift. The thermal diffusion of the dragged particles, as well as the application 
of an external colored noise, can markedly modify the direction and magnitude of 
the Stokes' drift velocity. Furthermore, it has been shown that such a stochastic 
Stokes' drift is equivalent to another paradigm of Brownian motors, the rocking 
ratchet |18j . In this scenario, the systematic motion acquired under the influence of 
an alternating zero-average stochastic force revolves essentially around the asymme- 
try of the nonlinear response in the presence of a steady asymmetric potential 20 . 
This mathematical equivalence between stochastic Stokes' drift and Brownian mo- 
tors hints at various potential applications of Stokes' drift, for instance, transport 
by capillary waves, storage of light in quantum wells, single-electron transport in 
one-dimensional channels, optical tweezing of colloidal particles, and diffusion of 
dislocations in solids, as noted in ^1 an d references therein. 

In this paper we show, with an analytically solvable model, that the character- 
istics of the one-dimensional stochastic Stokes' drift are quite complex, and that 
various interesting phenomena are induced, including enhancement of the deter- 
ministic drift and current reversal, when the particles are subjected to an additive 
colored noise, specifically, a dichotomous noise. The starting point is the following 
stochastic equation with an additive symmetric dichotomous perturbation: 

x = f(x- vt) +A£(t). (2) 

The dichotomous perturbation has an amplitude A and the stochastic variable £ (t) 
takes on the values ±1 with a transition rate k. It is appropriate to assume, on 
physical grounds, that the particle cannot move faster than the wave and thus 
f(y) < v for all y. The quantity of interest is the asymptotic drift velocity, (x) = 
limt_ >00 (x(f))/f, where the brackets indicate an average over the realizations of the 
dichotomous noise. Introducing a new variable y(t) — x(t) — vt, one can rewrite 
Eq. J5J as 

y = F(y) + M(t), (3) 

with a time-independent periodic force, F(y + L) = F(y) = f(y) — v < 0. Equa- 
tion © is a variant of the general stochastic differential equations solved in |§lllO| , 
with a dichotomous force F(y) ± A. 

The behavior of the system, and the corresponding solution of the associated 
master equation for the probability density, depend on whether or not there are 
unstable fixed points in the "±" dynamics. We present below the results for the 
two simplest cases, namely, one with no fixed points (in Sec. I2.l| l. and the other 
with two fixed points in the "+" dynamics and no fixed points in the "— " dynamics 
(in Sec.E3. 

2.1. Systems with no fixed points 

When A > max|F(y)| or < A < min|F(y)|, there are no fixed points. Then 
one obtains the following expressions (see Eqs. (7) and (8) in Ref. jTHj) for the 



asymptotic probability density, 
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and for the asymptotic mean velocity, 
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Here denotes the derivative with respect to y. The detailed analysis of 

these results, including the possibility of current reversal, were already discussed in 
Ref. QH. 

2.2. Systems with asymptotic fixed points 

If the amplitude of the noise lies in the intermediate regime min|_F(y)| < A < 
max then the "+" dynamics has at least one pair of fixed points in [0, L). 

We consider here the simplest case in which there is only one pair of fixed points, 
namely, a stable fixed point y\ with F'(yi) < 0, and an unstable fixed point y2 > yi 
with F'(y2) > 0. However, systems with several pairs of fixed points can be treated 
in the same way. According to the general discussion in Section III.B of Ref. |10|. 
the probability density in the interval [yi, y\ + L) is given by 
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Equation Q is continuous at the unstable fixed point y2, with 

lim P(y) = lim P(y) = _M f r +F '^\ r 
y/V2 yy > y\ V 2 yyl L 2A{k + F'(y 2 )} 



(7) 



Noting that F'(y2) > and P{y2) 0, (y) is necessarily zero or negative. This 
immediately suggests that the direction of the net drift (i) — (y) +v can be reversed 
by varying the amplitude A or the transition rate k of the dichotomous perturbation. 
The details of this current reversal phenomenon are discussed in a later section. At 
the stable fixed point y\ the probability density is continuous only for k/\F'{yi)\ > 1 
and its value is 



(8) 



For k/\F'(yi)\ < 1, P{y) is divergent at y\ but integrable. 
The corresponding mean velocity is given by 
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Note that if several waves are present, their contributions are not additive due to 
the highly nonlinear dependence of (x) on F. This seems to be a general feature of 
the stochastic Stokes' drift |HU16llT7| . contrary to the deterministic drift. 

2.3. The square wave 

The above results are valid for an arbitrary form of the wave F(y), provided that the 
obvious necessary differentiability and integrability conditions are fulfilled. How- 
ever, due to the multiple integrals in Eqs. JHJ and 0, further analytic investigation 
without a simple form of F(y) is difficult. Following our previous work |11U1(J| . we 
use a piecewise linear wave 
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(10) 



with < b < 1, and F(y + L) = F(y). Although the analytic treatment is possible, 
the general solutions are still rather complicated, and the results are relegated to 
|Appendix A| We focus here on the results in the limit of a square wave I — > 0, for 
which 

' -F- = -(1 - b)v for y e [0, L/2) 

(11) 

-F+ = -(l + b)v for ye [L/2, L) 



It should be noted that this limit is singular and must be handled with some care. 
From Eq. (@J one obtains the probability density for A > F+ or < A < F_ , 
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Fig 1. The asymptotic Stokes' drift velocity (normalized by its deterministic value vo = b 2 v) 
as a function of normalized transition rate /3 = 2kvL/A 2 for three different cases. The symbols 
are the results of numerical simulations. Circles: A = 0.4 < F— (no fixed point); squares: 
F- < A = 1.0 < -F+ (two fixed points); triangles: A = 1.6 > F+ (no fixed point). The solid lines 
indicate the results of theory. The values of the other parameters arc v = 1, b = 0.5, and L = 1. 
The direction of the drift velocity remains the same for A = 0.4 and A = 1.0, but is reversed for 
A = 1.6. 



Integrating Eq. ©, we find the probability for F_ < A < F + , 
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where 5±(x) are the half Dirac-<5 functions |21|. 

Correspondingly, the drift velocity for A > F + or < A < F_ is given by 
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where, recall, vq = b 2 v is the deterministic value of the Stokes' drift velocity for the 
square wave. We mention some limiting cases of interest, (i) The limit A — > or 
k — > oo leads to the deterministic Stokes' velocity, (i) = v . (ii) The drift velocity 
derived in Ref. |16| is recovered in the white noise limit A — > oo, k — ► oo, with a 
finite I? = A 2 /2fc. (iii) The quenched-noise limit k — » describes the system in 
which half of the particles, chosen at random, are subjected to a constant external 
forcing +^4 and the other half to forcing —A. When there are no fixed points, the 
mean velocity in this limit is 



<*> - v ° [y^A2 ) • (17) 

This expression clearly indicates the existence of the flux reversal suggested earlier, 
since the drift velocity is negative for sufficiently large amplitude A of the forcing, 
and positive for small A's. When there are fixed points in the "+" dynamics, 
Eq. Qlfijl leads to a mean velocity 

, v F 2 + F 2 - 2A 2 



In this limit, the current is reversed at A = y (F 2 + F 2 )/2 = \[\ + b 2 v. 

Figure n shows the drift velocity, Eqs. H15[) and (|lfci|) . as a function of the nor- 
malized transition rate, (3 = 2kvL/A 2 . The plots of the analytic results are in good 
agreement with the Monte Carlo simulation of the original stochastic differential 
equation @. When the transition rate is high (/3 > 1), the effect of the dichoto- 
mous noise is averaged out and the drift velocity approaches the deterministic limit 
Vq regardless of the noise amplitude. Interesting phenomena, namely, negative drift 
velocity and current reversal, take place only when the transition rate is sufficiently 
small, i.e., < 1. 

Figure|21shows the drift velocity, Eqs. (|15|l and (|l(j|> . in the three different regimes 
corresponding to low (A < PL) and high (A > F + ) noise amplitudes with no fixed 
point, and intermediate (PL < A < F + ) noise amplitudes with a pair of fixed 
points. The agreement between analytic results and simulation is again very good. 
At low noise amplitudes the drift velocity is positive and increases with increasing 
noise amplitude, i.e., noise enhances Stokes' drift. The noise-induced enhancement 
reaches its maximum at A = P_ and decreases above this noise amplitude. In 
general, noise-induced phenomena disappear when the noise becomes too large. 
However, the decrease at the intermediate noise is not the destruction of the Stokes' 
drift due to large fluctuations of the particle velocity. The sudden change with 
discontinuity in the first derivative of the drift velocity suggests the appearance of 
a bifurcation in the dynamics. When P_ < A < P + , the system has a stable fixed 
point in the "+" dynamics and some particles become stuck at the fixed point until 
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Fig 2. The asymptotic Stokes' drift velocity (normalized by its deterministic value vq = b 2 v) as a 
function of the noise amplitude A. The solid line represents theory and the circles show the results 
of numerical simulations. The values of the parameters are v = 1, b = 0.5, k = 1, and L = 1. 

the noise switches. The actual drift takes place only when the system is in the "-" 
dynamics, corresponding to a smaller or even negative value of x. The negative drift 
reaches its maximum at A = F + where the fixed points disappear. For A > F+, the 
"±" dynamics correspond to drift velocities in opposite directions, and therefore the 
net drift is reduced. Finally, as expected, Stokes' drift is asymptotically destroyed 
as the noise amplitude increases to infinity. 

Figure illustrates the probability densities Eqs. I|12|l. (|14|l for three different 
amplitudes of the dichotomous noise. For comparison, the numerical solution of 
the stochastic differential equation J3J is also shown. The agreement between the 
present analytic theory and the computer simulation is nearly perfect. Due to 
the discontinuity of the square wave F(y), Eq. I|1U|) . the probability densities are 
discontinuous at y — and y — L/2 for all cases. Moreover, for A — 1.0 the system 
has fixed points and the probability density becomes 5-singular at the stable fixed 
point. The probability densities are also asymmetric in all cases. For A = 0.4, the 
density is higher at y £ (0, L/2) than at y <E (L/2,L). On the other hand, the 
situation is the other way around for A = 1.6. This difference causes the current 
reversal. 

3. Hypersensitive Transport 

Generally and somewhat loosely speaking, the term "hypersensitive transport" 
refers to a highly nonlinear directed response of a nonequilibrium, noisy system to a 
small systematic external forcing. This phenomenon was discovered rather recently, 
and received a great deal of theoretical [§lll2| and experimental JSj attention. 
One of the simplest models exhibiting this novel phenomenon describes an over- 




Fig 3. The profile of the probability density P(y) for three values of the amplitude A of the 
dichotomous noise. The first (A < F-) and the last panel (A > F+) correspond to the regime 
with no fixed points, while the second panel (F— < A < F+) corresponds to the presence of two 
fixed points. Solid lines indicate the analytic results and shaded histograms show the results of 
numerical simulation. The values of the parameters are v = 1, 6 = 0.5, k = 1, and L = 1. In the 
simulation, the dynamics of 20000 particles is sampled at 100 different times. 
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Fig 4. The net potentials ^fV(x) — Fx. x\ and X2 represent unstable and stable fixed points, 
respectively. The substrate potential is defined as V(x) = — J x v(x')dx' . When the switching rate 
between these two potentials is sufficiently small, particles accumulate at a stable fixed point X2, 
bypassing the unstable fixed point x.\. When the potential switches, the particles are repelled by 
the unstable fixed point x\ toward a new stable fixed point. Even when the external force F is 
very small, the particles are able to avoid backward drift. 

damped particle whose dynamics switches dichotomically between a symmetric po- 
tential and its negative. When the symmetry of the system is slightly broken by a 
small directed external force, the system responds highly nonlinearly by exhibiting 
a "giant" systematic particle drift. Let us thus consider the following stochastic 
dynamics with a multiplicative noise 9 : 

x = F + £(t)v(x). (19) 

The dichotomous noise = ±1 has, as before, transition rate k. F > represents 
a constant external force, and v(x) is a given symmetric "substrate" force profile 
that is assumed to be periodic, v(x + L) — v(x). Some preliminary results on this 
model were reported in Ref. [§]. Here we present more (and more detailed) results 
than in our earlier work. 

3.1. No fixed points 

When the external force is sufficiently large [F 2 — v 2 (x) ^ for any x £ [0, L)], there 
is no fixed point in either dynamics, and one obtains the asymptotic probability 
density 



and, through the normalization of P(x), the mean asymptotic velocity 
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3.2. Asymptotic dynamics with fixed points 

For simplicity and without the loss of any relevant point of the method, we take 
v(x) to be a continuously decreasing function in [0, L/2] and symmetric about L/2, 
v(x + L/2) — —v(x). This implies that P(x + L/2) = P(x), so that we can limit 
our analysis to half a period. When the external force is weak, the particle moves 
alternately in the two "net potentials" represented schematically in Fig. 21 In this 
simple case the equation F 2 — v 2 (x) — has only two solutions in [0, L/2], x%, 
corresponding to an unstable fixed point in the "— " dynamics [F = v(x\), v'{x\) < 
0], and X2, a stable fixed point in the "+" dynamics [F = — v(x2), v'{x<2) < 0], with 
X2 > x\. These fixed points are the local extrema of the net potentials. 

According to the discussions in Pllltj , the physically acceptable solution in the 
interval \xi — L/2, X2] is given by 



{ > LF \ \F 2 - v 2 (x)\ 
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which extends by periodicity to the whole x-axis. At the unstable fixed point X\, 
the probability density is continuous and its value is given by 



lim P(x) = lim P(x) — — 

x\xi x/x; LF {1 
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Note that the stable fixed points X2—L/2 and X2 are located at the ends of the chosen 
interval and there is no other stable fixed point within this interval. Therefore, the 
probability density is continuous inside the interval. However, it can be singular at 
the stable fixed points. When the transition rate is sufficiently large [fc/|t/(a;2)| > 1], 
there is not enough time for the particles to reach and to accumulate at the stable 
fixed points, and thus P(x) is continuous X2 and X2 — L/2 as follows: 
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For low transition rates [fc/|u'(^2)| *S 1], the majority of particles become stuck near 
the stable fixed points for a long time, thus resulting in the divergence of P(x) at 
X2 and X2 — L/2. These divergences cause a highly nonlinear conductivity of the 
system, as discussed later. 



From the normalization of P(x), the average velocity is obtained as 
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Although exact, the above results are still too complicated to gain a clear phys- 
ical picture of the behavior of the system. We therefore turn to a particular shape 
of the velocity profile v(x) that simplifies the evaluation of the above functional 
expressions. 

3.3. Piecewise linear internal force 

We consider a piecewise linear "substrate" force 
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with I ^ L/A and the periodicity condition v(x + L) = v(x). It is convenient to 
introduce the following dimcnsionlcss variables: 



f = F/v , a = lk/v , r = 4l/L (0 < T < 1), 
and to work with the function 
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where \ = L/ (21) — 1 — x/7. 

In the absence of fixed points (/ > 1), one obtains the stationary probability 
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The corresponding mean velocity is given by 
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When < / < 1 the dynamics has two fixed points, an unstable one at x\ — 
L/2 — 1(1 + /) and a stable one at x 2 = L/2 — 1(1 — /). The situation is then 
more complicated but still analytically tractable. The probability density is written 
separately for three different regions: 
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This probability density is continuous at the unstable fixed point x\, 
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Fig 5. The mean asymptotic velocity as a function of the applied force / for different values of 
the parameter a. The other parameters are i>o = 1, T = 0.4, and L = 1. The symbols represent 
the results of numerical simulations. The solid lines are the result of theory and the dashed line 
indicates the linear response. 



At the stable fixed point x%, P(x) is continuous for a > 1, with 



lim P(x) = lim P(x) = — — - r tt—, rrr- 

x\ X2 y ' X /X2 v ' LF{l + l/[2(a- 1)]}' 



(34) 



but divergent and integrable for a ^ 1 . The average velocity as usual follows from 
the normalization of P(x), 



= <1- 



iaf 



1 + / 
1-/ 



2a 



A" 1 -1 



r(i + /)° 

2A(1 - f)° 



dt 



-1 
s 



dt 



t[T[L/2-l(t+l)]~T( Xl )] 



|*-/| 1+a |* + /| 1_Q 



[T( Xl )(l + A)-T(L/2)] (35) 



The integrals in Eqs. (|28p. jSJ , and l|35|) cannot be evaluated in closed analytic 
form except for a few specific values of a. As an illustration, in |Appendix B| we 
present closed analytic results for a — 1/2. 

3.4. Response Properties 

Equations l|3T)) and (|3*5|l are numerically evaluated and plotted in Fig.|3]as a function 
of the external force / for different values of a. The results are in good agreement 




Fig 6. The profile of the probability density P(x) for two values of the applied force /. The first 
corresponds to the regime with fixed points, and the second to the regime with no fixed points. 
Solid lines indicate the analytic results and shaded histograms show the results of numerical 
simulations. The values of the other parameters are vo = 1, ct = 1, V = 0.4, and L = 1. In the 
simulations, the dynamics of 20000 particles is sampled at 100 different times. 



with computer simulations. The drift velocity shows highly nonlinear responses to 
the external force depending on the transition rate of the dichotomous noise, as 
illustrated in Fig. [S] There are four asymptotic regimes of interest. 

(i) Linear Response Regime I: When the deterministic external force / is much 
larger than the fluctuating force (/ 3> 1), the effect of the fluctuating force is 
negligible and a linear response, (x)/vq — /, is expected for any transition rate. 
Indeed, all curves in Fig. [5] approach the linear response curve as / increases. 

(ii) Linear Response Regime II: When the transition rate is very high (a 3> 
1), the fluctuating force is averaged out and only the deterministic external force 
effectively drives the particles. Therefore, a linear response, (x)/vq = /, is again 
expected for a 3> f . Figure [S] shows that the drift velocity for a = 1 is already very 
close to this linear response limit. 

(hi) Adiabatic Regime: Since we are interested in nonlinear response, we focus 



on the cases with small a and small /. When the transition rate is sufficiently 
small (adiabatic limit), the system remains in one of the ±V(x) potentials for a 
long time. The particles eventually reach a minimum of the net potential and 
wait there until the potential flips. When the dichotomous noise switches its value, 
the particles suddenly find themselves near the maximum of the net potential as 
illustrated in Fig. They move down to a new potential minimum and wait again 
for the next potential flip. A typical time for them to escape from the region 
close to the maximum is given by r — —(£/vo)]nf. Therefore, this adiabatic 
regime is realized when the average time between switches, is much longer 

than r, which leads to 1 > / 3> exp(— 1/a). (Note that the adiabatic regime is 
not possible unless there is a fixed point and thus / < f .) Only q«1 can satisfy 
this condition. In this range of the external force, the mean velocity is simply 
half the spatial period of the "substrate" potential divided by the mean switching 
time, that is, (x) = Lk/2 = 2voa/T. The same result can be obtained directly 
from Eq. ll.''>5ll by taking the appropriate limit. This result is interesting in that 
the mean velocity does not depend on the external force F, nor the magnitude 
vq of the fluctuating force, and is inversely proportional to the correlation time 
of the noise. A higher transition rate reduces the waiting period at the potential 
minimum and thus increases the mean velocity. However, increasing k reduces the 
range of occurrence of this adiabatic regime. These results are in perfect agreement 
with the computer simulations (Fig. |SJ . The transport in the adiabatic regime has 
been investigated as hypersensitive response in Ref. |12| but without a full analytic 
solution. 

(iv) Hyper- Nonlinear Regime: For a small external force, / < exp(— 1/a), the 
particles manage to advance to the next potential minimum only when the dichoto- 
mous noise realizes exponentially rare cases where it keeps the same value for times 
much longer than the correlation time. Therefore, the mean velocity rapidly falls 
to zero as / decreases. Indeed, from Eq. 1|35JI ■ taking the limit / — ► 0, we find for 
the mean velocity 

(x)/v w parQn/) 2 ]- 1 . (36) 

This result is striking in that the susceptibility (or the conductance) of the 
system d(x)/df diverges at / = 0, indicating that the mean velocity is extremely 
sensitive to the external force. Another interesting feature is that the mean velocity 
is inversely proportional to the transition rate k, in contrast to the adiabatic regime 
where the velocity is directly proportional to the transition rate. Figure shows 
the mean velocity as a function of the applied force / for different values of a. Note 
the linear response regime I for all values of a; the linear response regime II for 
a = 1 (practically for all /, except the small-/ region); the adiabatic regime for 
a = 0.01 and (more restricted) for a = 0.1; and the hyper- nonlinear regime for 
all a. The nonmonotonic dependence of the mean velocity on a (e.g., ~ « in the 
adiabatic regime, ~ oT 1 in the hyper- nonlinear regime, and a-independent in the 
linear response regime) leads to the fact that the curves in Fig. [5] cross each other 
in the /-regions where there is a switch between these various regimes. 

In Fig. HO we present two typical profiles of the probability density P(x) for the 
two regimes with and without fixed points, with an abrupt change in shape between 
these regimes. 



4. Conclusions 



We have derived explicit results for the probability density and drift velocity in sys- 
tems that model Stokes' drift and hypersensitive response driven by dichotomous 
noise. We include the situation in which the asymptotic dynamics crosses unstable 
fixed points, a case where the standard approaches are not applicable in a straight- 
forward way HOES- We have presented analytic results for particular choices of 
potentials and parameters that elucidate the behavior in a way not entirely possible 
when only numerical results are available. 

Appendix A. Stokes' Drift Problem: Results for the Piecewise Linear 
Wave Eq. (fTU|) 

In the absence of fixed points one finds the following general piecewise expression 
for the probability density in four different regions: 
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(A.l) 

where we have used the dimensionless quantities a = lk/bv, and <j>\, (f>2 as defined in 
Eq. H13|l. Note that the system is invariant under the transformation, y — > (y — L/2) 
and b — > —6. Therefore, P(y) for y > L/2 can formally be obtained from the 
expression for P(y — L/2) by simply changing b to —b and {nota bene) a to —a. 
The integration constants C1-C4 assure that P(y) has no discontinuity in this case 
and is periodic. The mean asymptotic velocity (i) — v + (y) is obtained by imposing 
the normalization of P(y)- 

We turn now to the case of principal interest to us, F_ < A < F + , where the 
dynamics exhibits two fixed points, one stable, yi = L/2 — 21 + I [A — F-] /(bv), 
and one unstable, y2 = L — 21 + I [F + — A] / (bv). The piecewise expression for the 



probability density is now given by 
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The integration constants D\—D± are obtained by imposing the continuity of P{y) at 
y = L/2 - 21, L/2, and L - 21, and its periodicity. Note from Eq. (|A.2|I that P(y) is 
continuous at y 2 , ]im y ^ y2 P(y) = lim yNy2 P(y) = -({y)/L)[(2a + l)]/[2(a + 1)A]. 
Also, (y) < [from the positiveness of P(y 2 )]- At the stable fixed point, P(y) 
is either divergent, with a weak, integrable singularity whenever a ^ 1, or finite, 
lim^ yi P(y) = lim y ^ yi P(y) = -{{y)/L)[{2a - l)][2(a - 1)A] when a > 1. Fi- 
nally, the mean velocity is obtained by imposing the normalization condition for 

P{y)- 

The integrals that appear in the expressions ifA - !]) and fO|) of P(y) and, cor- 
respondingly, the integration constants C\— C4 and D±—Dn can be evaluated in a 
closed analytic form only for certain values of a. As an example, we present below 
these results for the particular case a = 1/2. 



When A > F + , the probability density reads 
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For < A < F_ , one obtains the probability density 
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and the corresponding asymptotic drift velocity 
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For F_ < A < F + , the probability density is expressed as 
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and the mean drift velocity as 
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where the integration constants D1-D4 are given by 
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Appendix B. The Problem of Hypersensitive Transport: Particular Cases 
for the Piecewise Linear Velocity Profile 

Again, explicit results can be obtained for the case a = 1/2. When there are no 
fixed points one obtains for the probability density: 
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where x = 2T — 1 — a;/Z. The mean velocity is 



(B.12) 
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When there are two fixed points, P{x) has an integrable power-law singularity 



at the stable fixed point: 
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For the mean velocity we find 



(x) 



= <U + r 



1 4 + 2 2 + 2 P ^ T(l - p) 



+ 



1 — exp I — 



2/(1 -r) 
r(i-/2) 



ryr+7 



/In 




/ 



i + ^ cxP r 2/(1 - r) ' 




(B.16) 



- (l + 2/)ji±|exp 



Acknowledgments 

This work was partially supported by the Swiss National Science Foundation (LB.) 
and by the National Science Foundation under grant Nos. PHY-0354937. 



References 

[1] F. Moss and P. V. E. McClintock, eds., Noise in Nonlinear Dynamical Systems, 
Cambridge University Press, Cambridge (1989); K. Lindenberg, B. J. West, and G. 
P. Tsironis, Bistable systems driven by colored noise, Rev. Solid State Sci. 3 (1989) 
143; K. Lindenberg and B. J. West, The Nonequilibrium Statistical Mechanics of 
Open and Closed Systems, VCH Publishers, New York (1990); P. Hanggi and P. 
Jung, Colored noise in dynamical systems, Advances in Chemical Physics LXXXIX 
(1995),239. 

[2] J. M. Sancho, External dichotomous noise: The problem of the mean-first-passage 
time, Phys. Rev. A 31 (1985) 3523; J. Masoliver, K. Lindenberg, and B. J. West, 
First-passage times for non-Markovian processes, Phys. Rev. A 33 (1986) 2177; J. 
Masoliver, K. Lindenberg, and B. J. West, First-passage times for non-Markovian 
processes: Correlated impacts on a free process, Phys. Rev. A 34 (1986) 1481; J. 
Masoliver, K. Lindenberg, and B. J. West, First-passage times for non-Markovian 
processes: Correlated impacts on bound processes, Phys. Rev. A (1986) 2351; M. A. 
Rodriguez and L. Pesquera, First-passage times for non-Markovian processes driven 
by dichotomic Markov noise, Phys. Rev. A 34 (1986) 4532; C. R. Doering, Comment 
on first-passage times for processes driven by dichotomous fluctuations, Phys. Rev. A 
35 (1987) 3166; V. Balakrishnan, C. Van den Broeck, and P. Hanggi, First-passage 
times of non-Markovian processes: The case of a reflecting boundary, Phys. Rev. A 
38 (1988) 4213; U. Behn and K. Schiele, Stratonovich model driven by dichotomous 
noise: Mean first passage time, Z. Phys. B 77 (1989) 485; M. Kus, E. Wajnryb 
and K. Wodkiewicz, Mean first-passage time in the presence of colored noise: A 
random-telegraph-signal approach, Phys. Rev. A 43 (1991) 4167; J. Olarrea, J. M. R. 



Parrondo, and F. J. de la Rubia, Escape statistics for systems driven by dichotomous 
noise. I. General theory, J. Stat. Phys. 79 (1995) 683; J. Olarrea, J. M. R. Parrondo, 
and F. J. de la Rubia, Escape statistics for systems driven by dichotomous noise. II. 
The imperfect pitchfork bifurcation as a case study, J. Stat. Phys. 79 (1995) 669; 
C. Doering, W. Horsthemke, and J. Riordan, Nonequilibmum fluctuation-induced 
transport, Phys. Rev. Lett. 72 (1994) 2984. 
[3] S. Mangioni, R. R. Deza, H. S. Wio, and Toral, Disordering effects of color in 
nonequilibrium phase transitions induced by multiplicative noise, Phys. Rev. Lett. 
79 (1997) 2389. 

[4] S. E. Mangioni, R. R. Deza, R. Toral, and H. S. Wio, Nonequilibrium phase tran- 
sitions induced by multiplicative noise: Effects of self- correlation, Phys. Rev. E 61 
(2000) 223. 

[5] Y. Kim and W. Sung, Does stochastic resonance occur in periodic potentials?, Phys. 

Rev. £ 57 (1998) 6237(R). 
[6] R. Rozenfeld, A. Neiman, and L. Schimansky-Geier, Stochastic resonance enhanced 

by dichotomic noise in a bistable system, Phys. Rev. E 62 (2000) 3031 (R). 
[7] R. Rozenfeld, J. A. Freund, A. Neiman, and L. Schimansky-Geier, Noise-induced 

phase synchronization enhanced by dichotomic noise, Phys. Rev. EGA (2001) 051107. 
[8] N. G. Van Kampen, Stochastic Processes in Physics and Chemistry, North-Holland, 

Amsterdam (1992). 

[9] I. Bena, C. Van den Broeck, R. Kawai, and K. Lindenberg, Nonlinear response with 

dichotomous noise, Phys. Rev. E 66 (2002) 045603(R). 
[10] I. Bena, C. Van den Broeck, R. Kawai, and K. Lindenberg, Drift by dichotomous 

Markov noise, Phys. Rev. E 68 (2003) 041111. 
[11] I. Bena, M. Copelli, and C. Van den Broeck, Stokes' drift: A rocking ratchet, J. Stat. 

Phys. 101 (2000) 415. 

[12] M. B. Tarlie and R. D. Astumian, Optimal modulation of a Brownian ratchet and 
enhanced sensitivity to a weak external force, Proc. Natl. Acad. Sci. USA 95 (1998), 
2039; V. Berdichevsky and M. Gitterman, Stochastic resonance and ratchets - new 
manifestations, Physica A 249 (1998) 88; S. L. Ginzburg and M. A. Pustovoit, Hy- 
persensitive transport in a phase model with multiplicative stimulus, Phys. Lett. A 
291 (2001) 77; ibid, Noise-induced hypersensitivity to small time- dependent signals, 
Phys. Rev. Lett. 80 (1998) 4840; ibid, Hypersensitivity of a nonlinear system with 
multiplicative colored noise to an external periodic system, J. Exp. Ph. Phys. 89 

(1999) 801; ibid., Stochastic resonance, on-off intermittency and "physical mathe- 
matics" Europhys. Lett. 45 (1999), 540; O. V. Gerashchenko, S. L. Ginzburg, and 
M. A. Pustovoit, Hypersensitivity to small signals in a stochastic system with mul- 
tiplicative colored noise, Eur. Phys. J. B 19 (2001) 101; R. Mankin, A. Haljas, R. 
Tammelo, and D. Martila, Mechanism of hypersensitive transport in tilted sharp 
ratchets, Phys. Rev. E 68 (2003), 011105. 

[13] O. V. Gerashchenko, S. L. Ginzburg, and M. A. Pustovoit, Multiplicative-noise- 
induced amplification of weak signals in a system with on-off intermittency, J. Exp. 
Theor. Phys. Lett. 67, (1998), 997; ibid., Experimental observation of noise-induced 
sensitivity to small signals in a system with on-off intermittency, Eur. J. Phys. B 15 

(2000) , 335; S. L. Ginzburg and M. A. Pustovoit, Stochastic resonance in two-state 
model of membrane channel with comparable opening and closing rates, Phys. Rev. E 
66 (2002), 021107; ibid., Bursting dynamics of a model neuron induced by intrinsic 
channel noise, Fluct. Noise Lett. 3 (2003) L265. 

[14] A few recent papers from a rich bibliography on the subject include: A. M. Bratseth, 
On the estimation of transport characteristics of atmospheric data sets, Tellus 50A 



(1998) 451; S. M. Cox, Onset of Langmuir circulation when shear flow and Stokes 
drift are not parallel, Fluid Dyn. Res. 19 (1997) 149; H. C. Graber, B. K. Haus, 
R. D. Chapman and L. K. Shay, HF Radar comparisons with moored estimates of 
current speed and direction: Expected differences and implications, J. Geophys. Res. 
- Oceans 102 (1997) 18749; J. C. McWilliams, P. P. Sullivan and Chin-Hoh Moeng, 
Langmuir turbulence in the ocean, J. Fluid Mech. 334 (1997) 1; M. Nordsveen and 
A. F. Bertelsen, Wave induced secondary motions in stratified duct flow, Int. J. 
Multiphase Flow 23 (1997) 503; V. Polonichko, Generation of Langmuir circulation 
for nonaligned wind stress and the Stokes drift, J. Geophys. Res. 102 (1997) 15773; 
S. A. Thorpe, Interactions between internal waves and boundary layer vortices, J. 
Phys. Oceanogr. 27 (1997) 62; ibid., On the interactions of internal waves reflecting 
from slopes, (1997) 2072. 

[15] See e.g. C. P. Lee, The balancing of thermocapillary flow in a floating zone by ripple- 
driven streaming, Phys. of Fluids 10 (1998) 2765 and references therein. 

[16] C. Van den Broeck, Stokes' drift: An exact result, Europhys. Lett. 46 (1999) 1. 

[17] K. Herterich and K. Hasselmann, The horizontal diffusion of tracers by surface waves, 
J. Phys. Oceanogr. 12 (1982) 704; O. N. Mesquita, S. Kane and J. P. Gollub, Trans- 
port by capillary waves: Fluctuating Stokes drift, Phys. Rev. A 45 (1992) 3700; K. 
M. Janson and G. D. Lythe, Stochastic Stokes drift, Phys. Rev. Lett. 81 (1998) 3136; 
ibid., Stochastic calculus: Application to dynamic bifurcations and threshold cross- 
ings, J. Stat. Phys. 90 (1998) 227; R. Landauer and M. Biittiker, Drift and diffusion 
in reversible computation, Physica Scripta T9 (1985) 155. 

[18] For an overview, see P. Reimann, Brownian motors: noisy transport far from equi- 
librium, Phys. Rep. 361 (2002) 57 and references therein. 

[19] F. Marchesoni and M. Borromeo, Stokes drift of linear defects, Phys. Rev. B 65 
(2002) 184101. 

[20] C.R. Doering, W. Horsthemke and J. Riordan, Nonequiltbrium fluctuation-induced 
transport, Phys. Rev. Lett. 72 (1994) 2984; K. W. Kehr, K. Mussawisade, T. Wich- 
mann and W. Dieterich, Rectification by hopping motion through nonsymmetric po- 
tentials with strong bias, Phys. Rev. E 56 (1997) R2351. 

[21] The half delta functions are defined as 



lim 

e\,0 



J — oo 
/•+oo 



/ e 5-(x)f(x)dx = (1/2) \im f(x) = (l/2)/(-0) 





with 5(x) = 6- (x) + S+ (x) . 



